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ABSTRACT 



Context. Understanding massive star formation requires comprehensive knowledge about the initial conditions of this 
process. The cradles of massive stars are believed to be located in dense and massive molecular clumps. 
Aims. In this study, we present an unbiased sample of the earliest stages of massive star formation across 20 deg 2 of 
the sky. 

Methods. Within the region 10° < 1 < 20° and |bj < 1°, we search the ATLASGAL survey at 870 ^m for dense 
gas condensations. These clumps are carefully examined for indications of ongoing star formation using YSOs from 
the GLIMPSE source catalog as well as sources in the 24 /im MIPSGAL images, to search for starless clumps. We 
calculate the column densities as well as the kinematic distances and masses for sources where the v; sr is known from 
spectroscopic observations. 

Results. Within the given region, we identify 210 starless clumps* with peak column densities > 1 x 10 23 cm" 2 . In 
particular, we identify potential starless clumps on the other side of the Galaxy. The sizes of the clumps range between 
0.1 pc and 3 pc with masses between a few tens of Mq up to several ten thousands of Mq. Most of them may form 
massive stars, but in the 20 deg 2 area we only find 14 regions massive enough to form stars more massive than 20 
M© and 3 regions with the potential to form stars more massive than 40 Mq. The slope of the high-mass tail of the 
clump mass function for clumps on the near side of the Galaxy is a = 2.2 and, therefore, S alp eter- like. We estimate 
the lifetime of the most massive starless clumps to be (6 ± 5) x 10 4 yr. 

Conclusions. The sample offers a uniform selection of starless clumps. In the large area surveyed, we only find a few 
potential precursors of stars in the excess of 40 Mq . It appears that the lifetime of these clumps is somewhat shorter 
than their free-fall times, although both values agree within the errors. In addition, these are ideal objects for detailed 
studies and follow-up observations. 

Key words. Surveys, Stars: formation, massive, ISM: clouds, Infrared: stars, Submillimeter: ISM 
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1. Introduction 

The understanding of high-mass star formation has made 
tremendous progress within the last decade (cf. reviews by 
|Zinnecker fc Yorke|200~7] |Beuther et al.|2007[ ) . Nevertheless , 
there is still no consistent scenario describing how massive 
stars form, nor is the impact of massive stars on their en- 
vironment fully understood. They undoubtedly play a cru- 
cial role throughout their whole lifetime on the physical and 
chemical structure of the interstellar medium, the evolution 
of stellar clusters, and even of entire galaxies. Therefore, 
massive star formation needs to be understood if we wish 
to make progress in our understanding of galaxy evolution. 

Indepen dent of the actual high-mass star formation sce- 
nario (e.g. |Keto[ |2003[ |McKee fc Tan] |2003[ |Bonnell fc 



The catalog (Table \3\ is available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 



via http://cdsweb.u-strasbg.fr/cgi-bin/qcat '/J/ A+A/ 



Bate 2005 Commergon et al. 2011), there is agreement 
that the most massive stars form in clusters. Therefore, it 
is probable that we can detect an initial stage of a high- 
mass gas clump without any star formation signatures. 



The discovery of infrared-dark clouds (IRDCs) with 
ISO, MSX, and Spitzer data provided an interesting sample 
of objects with which to ch aracterize the earliest stages of 



mass ive star formation (e.g. Perault et al. 1996 



1998). The systematic study of IRDCs provided potential 



Carey et al. 



precursors of massive stars and allowed the chara cteriza- 
tion of their physical and chemical parameters (e.g. Simon 
et^[2006l |Peretto fc Fuller ||2Tj09l fVasyunina et al | |2009, 
20l0j . Although |Peretto fc Fullerl ( |2009[ ) reported that more 
than 30 % of the IRDCs have no IR counterparts at 24 /im, 
follow-up studie s of IRDCs often revea l ed signs of ongoing 



star formation (Beuther fc Sridharan[ 2007 Cyganowski 



et al.[|2008[|Vasyunina et al.[|2010 



iranl 
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All IRDC studies are biased by the variation of the 
background, foreground confusion, and extinction caused 
by variations in the dust properties. Longward of 200 /im, 
the thermal emission from dust grains in IRDCs is optically 
thin and can be measured at mm and sub-mm wavelengths 



(Hildcbrand 1983). This can be used to obtain extinction- 



or 10 mJy in either the 3.8, 4.5, 5.8, or 8.0 /j,m band, re- 
spectively, for a source to be taken into account. For the re- 
gion given, we investigated more than 5.6 million GLIMPSE 
sources. The IRAC/SPITZER pointing accuracy is better 
than 1" and the pixel resolution is 0.6". As an additional 
tracer of ongoi ng star formation, w e used the MIPS GAL 



independent mass measurements of the cold gas inside these 24 /xm survey ( Carey et al. 2009 1 . The rms noise of the 
objects. 

Until recently, there were no available systematic sur- 
veys of the Galactic plane directly tracing the cold dust 
associated with molecular clumps. Now, the Bolocam 



Galactic Plane Survey (BGPS, jAguirre et al.||2011[ ) at 1.1 
mm, as well as the APEX Telescope Large Area Survey 
of the GALaxy (ATLASGAL, |Schuller et aL][2009| at 870 
fim, offer surveys of the Galactic plane's cold dust. Only 
ATLASGAL covers the full range 1 = -60° to 60° of the 
inner Galactic plane at 19" resolution. 

In this paper, we present a compilation of clumps of high 
column density, located in a region of Galactic longitude 10° 
< 1 < 20° and latitude -1° < b < 1°, showing no signs of star 
formation. To confirm their starless nature, we carefully 



examined each clump for GLIMPSE (Benjamin et al. 2003 1 



and MIPS GAL 24 /<m ( |Carey et al.||2009| ) sources, eTthe: 
of which would indicate that star formation had already 
started. The column density threshold we imposed on our 
survey is 1 x 10 23 cm -2 . 

After a short description of the surveys we employed 
(Sec. |2|, we describe the individual steps of the classifica- 
tion and its limitations in Sec. [3} In Sec. [4] we present both 
the results of the classification and the direct clump prop- 
erties and discuss the clump column densities. Using the 
ammonia velocities given in Wienen et al. (submitted), we 
derived distances to ~ 71 % of the clumps (Sec. [5]). Sec. [6] 
presents the clump masses and the clump mass function. 
Next we estimated the lifetimes of starless clumps (Sec. u\. 
In Sec. [8] we discu ss our results and compare them to other 
surveys (Sec. 8.3). Our conclusions, in Sec. |9j summarizes 
the main results of this work and provides an outlook to 
future work. 



2. Employed data 

All data for this study were taken from large surveys, 
most of them publicly available. Clumps were identi- 
fied by searching for continuum peaks at 870 /xm in the 



ATLASGAL survey (Schuller et al. 2009) and then clas 



sificd using both the GLIMPSE point source catalog and 
MIPSGAL 24 /xm images. 

In contrast to most other searches for massive prestel- 
lar clumps of high column density using ext i nction maps 



([Simon et al.[ |2006| |Peretto fc Fuller) |2009| |Kainulainen 



et al.[ 201ip , we used emission at 870 /xm as a tracer of cold 



dust. The AP EX telescope large ar ea survey of the Galaxy 
(ATLASGAL, |Schuller et al.||2009[ ) is a systematic survey 
of the Galac tic plane at 870 /mi with LABOCA (Siringo 
et al.[ |2009|). Its beam size is 19.2", the pixel size in the 



maps is 6" , and the average rms noise of the selected maps 
is below 50 mJy. To obtain a statistically meaningful sam- 
ple, we covered 20 deg 2 on the sky, the region of Galactic 
longitude 10° < 1 < 20° and Galactic latitude |b| < 1°. 

As described in Sec. |3.2[ we extracted young stellar ob- 
jects f rom t he GLIMPSE I Spring '07 catalog (Benjamin 



et al. 2003). Among other criteria, the GLIMPSE source 



MIPSGAL images is ~ 0.67mJy and its resolution is 6". 

3. Classification 

The naming of clumps in the literature refers to various 
physical objects. Research groups working on high- and 
low-mass star formation have different naming schemes for 
the objects named clumps and cores, including t he sub- 
categories starless and prestellar ( |Enoch et al. 2008 ) . An of- 
ten used nomenclature tries to bind things to physical prop- 
erties, denoting gravitation ally bound objects "cores" an d 
unbound objects "clumps" (Chabrier & Hennebelle 2010). 
In this paper, we use the term clumps for all emission peaks 



revea led by the CLUMPFIND algorithm (Williams et al 
19941). Typically, these are massive and large enough to 



form massive clusters. As shown in Fig. [5j typical sizes de- 
rived from the effective radii of these clumps range from 18" 
to 70". These can either be bound or unbound systems, but 
are assumed to be coherent in lbv space (Galactic longitude 
and latitud e, and radial velocity, thus distance, Williams 



et al. 2000). In the remainder of this paper, clumps are 
called starless if they host no mid-IR tracers of ongoing 
star formation. Nevertheless, as mentioned in Sec. [I] and 
Sec. |3.3.3| many IRDCs not hosting 24 /im sources have 
been shown to host star formation using other tracers such 
as SiO emission. In this study, we cannot completely rule 
out the presence of star formation, but only present star- 
less clump candidat es. In this context, the MALT90 survey 
(Foster et al. 2011) will improve future classifications. 



3.1. Clump extraction 

To identify starless clumps , we first emp l oyed the 
CLUMPFIND algorithm by |Williams efH] ([1994]) to 
search for dust condensations. It has been reported that 
CLUMPFIND is less reliable in very crowded regions 



(Kainulainen et al. 2009) and that the ext racted clump 
param eters strongly depend on the distance ( |Smith et al. 
2008). The second point is unavoidable in observed data 
and is discussed further in Sec. 16.21 Nevertheless, we are 



catalog requires a minimum flux of 0.6 mJy, 0.4 mJy, 2 mJy, 



interested in column-density peak positions and the associ- 
ated fluxes/masses , whic h CLUMPFIND can extract reli- 
ably. Pineda et al. (2009) demonstrated that the exponent 
of the derived mass function is not very sensitive to the 
chosen step size. 

We set the lowest detection level to 6<r, or 0.3 Jy. The 
additional thresholds, 0.3, 0.4, 0.5, 0.7, 0.9, 1.3, 1.8, 2.5, 4, 
and 7 Jy, were chosen (1) to account for the degree of varia- 
tion relative to the actual emission level and (2) to trace the 
structures recognized by observers. The use of non-constant 
intervals for the various emission levels in CLUMPFIND 
prevents bright clumps from being artificially sub-divided 
because of brightness changes that are very small relative 
to the flux level of the clump. 

To test the robustness of the chosen thresholds, we 
compared the integrated clump fluxes of our clump ex- 
traction to classical 3 a spaced thresholds as proposed in 
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11.2 11.1 11 0.011 0.01050.009 0.0085 0.008 

Galactic longitude projection onto longitude axis 



Fig. 2. The figures show the clump definition as used in this paper (bottom row), compared to 'classical' 3cr spacing 
contour levels (top row). While the left most column shows the ATLASGAL image of Gil. 11 with the clump definitions 
(red), the other columns show two profiles along the lines shown in the left panel ('cut A' and 'cut B'). The red lines 
indicate the clump borders. 



Kainulainen et al. ( |2009[ ). The two right columns of Fig. 
2] compare the flux distributions along the lines plotted in 
the left column. The upper row shows our clump definition, 
the lower row shows the clumps of a pure 3<r spacing of the 
contour levels. While some regions do not differ at all, e.g. 
as shown by the right-most column of Fig. [2j the additional 
contours in the evenly spaced situation can subdivide large 
clumps into a number of smaller ones, shown by the middle 
column of Fig. [2] 

While flux level spacings that are not bound to some 
objective criteria introduce subjectivity, Fig. [2] shows that 
the chosen levels trace structures that we call clumps. In 
addition, we compared the fluxes of our clumps to the 
fluxes of the 3 a extraction and could not find an excess of 
bright clumps, as might have been expected. Furthermore, 
we compared our CLUMPFIND sources to sources found 
by Contrera s (priv. communicatio n) using SExtractor as 
described in Schuller et al.| ( |2009| . For fluxes above our 
thresholds, almost all sources identified by Contreras had a 
counterpart within our clumps with matching peak fluxes. 
In addition to this, the comparison shows that we identified 
smaller fragments of lower mass, and the integrated fluxes 
in our catalog are lower than the integrated fluxes of the 
corresponding SExtractor clumps. This assures us that we 
do not produce sources with artificially high fluxes. 

As we aim to study sites of high-mass star formation, 
clumps with a peak flux of less than 0.5 Jy, correspond- 
ing to a column density of 1 x 10 23 cm -2 , are ignored in 
the following. (For a further discu ssio n of the derivation of 
column densities, we refer to Sec. 4.2 ) 



To motivate this threshold, we used the Orion neb- 
ula cluster (ONC). Its stellar m ass is about 1800 M Q 

To be consistent with 
we assumed a star for- 



1998). 

m 



(Hillenbrand & Hartmann 
estimates carried out in Sec 

mation efficiency of 30% and, therefore, estimated the ini- 
tial gas mass of the ONC to be 6000 M Q . As the cluster 
has dispersed over its lifetime, we set its initial radius to 
the typical radius w e fou nd for our most massive clumps, 
hence 0.7 pc (cf. Sec. 7.1 ). With the assumptions made and 
a spherically symmetric mass distribution, the initial peak 
column density in the ONC has been 1.8 x 10 23 cm -2 , or 
0.6 g cm" 3 . This agrees with the theoretical values found 
by Krumholz & McKee (2008). To avoid fragmentation in 
high-mass star formation, they require column densities of 
1 g cm" 2 , or 3xl0 23 cm" 2 . 



3.2. Identification of starless clumps 

Although it is unclear whether high- or low-mass stars form 
first, starless clumps should not host young stellar objects 
(class I sources, YSOs). To identify clumps hosting YSOs, 
we searched the GLIMPSE source catalog for stars with 
colors similar to known YSOs and compared those to our 
clumps. To do so , we fo llowed the classification given by 
Gutermuth et al.|(2008[). W e used additional color criteria, 
given in |Gutermuth et al" (2008) as well, to reject contam- 
inating extragalactic sources and AGNs 'that masquerade 
as bona fide YSOs'. Afterwards we selected objects obeying 
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starless clumps 



GLIMPSE catalog 

YSO / class I 



MIPSAL images 

warm dust 



clumps 

peak column density > 10 23 cm" 



visual inspection 

of starless 
clumps candidates 



starless clump candidates 




3.3. Limitations and observational biases 

Although the visual verification of the classification ensures 
a maximum reliability, technical limitations of the data sets 
impose various biases. To point out the limitations of this 
study, next we carefully discuss the biases. 



3.3.1. ATLASGAL and clump finding limitations 

The spatial limitations of ATLASGAL vary with the dis- 
tance and are discussed in detail in Sec. 16. 21 The flux thresh- 
old for the clump extraction was chosen to be ~ 6<r or 0.3 
Jy, and the higher thresholds were chosen to trace clearly 
recognizable structures. These threshold spacings are larger 
than the estimated rms. Contours in steps of the noise level 
are less biased and would generate more substructure, hence 
clumps. However, noise wo uld, more likely, g enerate artifi- 
cial clumps as studied by Reid et al. (20101. The chosen 
peak flux threshold of 0.5 Jy/beam corresponds to 1 x 
10 23 cm~ 2 . In the context of massive star formation, this 
provides a rough lower limit for p otential regions of mas- 
sive star formation (see Sec. 3.1). Sources for which the 



integrated flux is less than its peak flux are considered as 
artificial and 28 out of 929 sources were rejected. 



search peaks 
using clumpfind 




Fig. 1. Schematic visualization of classification. 



the following IRAC criteria: 

[4.5] - [5.8] > 1.0 OR 
([4.5] - [5.8] > 0.7 AND [3.6] 



[4.5] > 0.7). 



In addition, we required a source to be detected at 8 /jm. 
The identified YSOs were directly compared to the clumps 
and their extensions according to CLUMPFIND using the 
CLUMPFIND maps produced. If a YSO is located on a 
clump (in projection), the clump was considered as star 
forming and is ignored in the following. Nevertheless, at 
the onset of star formation, sources may be too cold to 
be detectable in the GLIMPSE bands, but show weak 24 
(im emission. Unfortunately, no MIPSGAL 24 /im point 
source catalog has been pub lished. Therefore, we used the 
STARFINDER algorithm by |Diolaiti etaL] ( |2000| to search 
the 24/im MIPSGAL images for point sources. To avoid 
mis-identifications we only extracted stars with a detection 
better than 7a. Again, clumps with a 24 /Ltm source were 
assumed to host stellar activity. In a last step, all remaining 
clumps were classified by visual inspection. Here the main 
focus was on 24 /im objects that had not been identified by 
STARFINDER. A schematic summary of the classification 
is given in Fig. [l] Parts of M17, in which MIPSGAL is 
saturated due to extended emission, were omitted as well as 
a few additional regions. Exact positions of omitted regions 
are listed in the Appendix, Table [B~T| 



3.3.2. GLIMPSE catalog limitations 

The detection thresholds of GLIMPSE and MIPSGAL as 
well as the point source extraction from the GLIMPSE 
catalog and the color-color criteria for young sources de- 
scribed in Sec.|3.2|have a major impact on the classification. 
Gutermuth et al. (20081 included criteria to reject contam- 
inating extragalactic sources, but AGB stars have similar 
colors to YSOs and obey the color-color criteria Gutermuth 
used to identify YSOs. Their co ntribution to the list o f 
YSOs may be as high as - 30% ( |Robitaille et al.[ |2008[ ), 
rejecting potential starless clumps. Nevertheless, their like- 
lihood of being projected onto a clump is significantly lower. 
Since we expect embedded YSOs to have detectable 24 fim 
flux for which we will check again, the given sensitivity 
limits of the GLIMPSE source catalog do not influence the 
population of starless clumps. Additional sources in the list 
of YSOs as well as chance alignments could lead to an ar- 
tificial rejection of clumps, but will not produce artificial 
starless clumps. 



3.3.3. MIPSGAL 24 /im limitations 

The situation is different for the 24 /im MIPSGAL im- 
ages. Here the sensitivity limit is the key parameter and 
sources hidden in the rms can lead to misidentifications 
of starless clumps. The brightness of the faintest sources 
still detectable varies over the images with respect to their 
surroundings, but for the visual inspection method we es- 
timated it to be ~ l mJy. This is slight ly smaller than th e 
2 mJy level for a 3tr detection given in Carey et al. (2009). 



For sources hidden in the dust, one may assume that 
all flux gets re-emitted by the dust producing a black body 
spectrum. From this, one can estimate the integrated lumi- 
nosity of the internal source. Since the faint 24 /jm sources 
in question are not detected at GLIMPSE wavelengths, we 
used both the GLIMPSE and MIPSGAL detection limits 
to construct a SED of three data points at 3.8 /zm, 8.0 /jm, 
and 24 ^m with 0.6 mJy, 10.0 mJy, and 2.0 mJy, respec- 



4 



J. Tackenberg: Starless clumps in ATLASGAL 



-0.36 



-0.38 



-0.4 
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Galactic longitjde / [ degrees ] 



10.55 



Fig. 3. Clump 67 overplotted on a MIPSGAL 24 /im im- 
age. The clump definition is in red and its peak position 
is mark ed by a red asterisk. Ov erplotted in green are HI I 



regions (Purcell & Hoare 



found in GLIMPSE by Robitaille et aLlpOOl. 



2010!), and in blue Red Sources 



tively. Using near- and mid-IR dust opacities from |Draine| 
& Lee (1984), we fitted a black body spectrum to the SED, 
with an integrated luminosity of 1.1 L at 3 kpc, or 27 L 
at 15 kpc. These luminositie s correspond to ma i n sequence 
stars of 1.1 M Q or 2.1 M Q (|Siess et al.j |2000[). |Krumholz 



et al. (20071 found that accretion luminosities in massive 



star formation reach several hundred solar luminosities very 
early on, excluding that a massive collapsing core could be 
hidden in the dust. Furthermore, there have been observa- 
tions of objects tha t may form high-mass stars, but are not 



yet that luminous (Beuthe r fc Steinacker 2007 Bontemps 
et al. 2010a Motte et al. 2010[ Ragan et al., in prep 



Nevertheless^ their luminosities are still higher (on the or- 
der of several 10 L ) than our detection limit for 24 /im 
fluxes on the far side of the Galaxy. Therefore, only low- 
mass objects can be hidden. 

Observations using the Herschel satellite have shown 
that some 24 /im dark regions, hence starless clumps al- 



(Beuther et al-l 2010 


Wilcock 


l Henning et al. 


(|20] 


M, these 



(2007) and Russeil et al. (2010) find MSX and 



Motte et al [ 

24 /im dark core s, driving SiO outflo ws. Although |Motte| 
et al. (20071 and Russeil et al. (2010) are less sensitive at 
24 /im, future studies will need to disentangle this situation. 



3.4. Verification of classification by comparison to other 
tracers and studies 



To test the classification 
and catalogs i s helpful. 
Similar to 



a comparison with other tracers 



(2008), Robitaille 



et al. 
tEe 



Gutermuth et al. 
(20081) identified Intrinsically Red Sources (R08 m 
following) by applying color-color criteria to the Spitzer 
GLIMPSE catalog. Both sets of color criteria differ and, 
in addition to a large number of common sources, both 
catalogs also identify different sources. We take these dif- 



ferent identifications as statistical variations that set the 
approximate uncertainties in the different catalogs. If we 
now compare the population of starless clumps to the Red 
Sources given in R08, this gives us a feeling for the classi- 
fication statistics. As it turns out, only two clumps that we 
identified as starless have a SPITZER Red Source. With 
knowledge of its position, we have been able to associate 
the R08 source in clump 67 with a peak in the 24 /im im- 
age (see Fig. [3]). The other Red Source is supposed to be 
in clump 1216, at a ridge of bright continuous emission. 
This hampers the identification and we cannot identify a 
24 /im counter part. Therefore, it is unclear to us whether 
this source is still very young. 

Another test of our classification is to check the clumps 
for additional tracers of star formation. Hn regions are a 
well-accepted tracer of ongoing massive-star formation and 
several surveys have searched the Galactic plane systema ti- 
cally. In this context, CORNISH ( |Purcell fc Hoare] |2010[ ), a 
Galactic plane survey at 5 GHz with the VLA in B config- 
uration, identified more than 600 Hn regions in our region 
of study. With their high spatial resolution, matches can be 
made unambiguously. We found Hn regions on only three 
of our clumps. In clump 67 (shown in Fig. [3| and clump 
87, no 24 /jm source is in the vicinity of the cm emission 
peak. This suggests that star formation is already taking 
place, but so embedded that (almost) no light can escape. 
However, clump 505 has a very bright 24 /im source at the 
edge of the clump. This might power the Hn region, which 
is offset by 11" towards the emission peak of clump 505. 

A comparison to t he Green Bank Te lescope Hn Region 
Survey (GBT HRDS, |Bania et al.[f20Tol and the Red MSX 
Sourc e Survey (RMS, ^Hoare et a ' 



2004 Mottram et al. 



2011) did not discover any matches 

In summary, since only three incorrect classifications 
have been found, all tests confirm our classification and 
establish its credibility. For consistency reasons, we flagged 
clump 67 as star forming, but kept the other two sources 
in our sample. 

Furthermore, both pointed HERSCHEL observations 
(e.g. EPOS, Krause et al., in prep., Ragan et al., in 
prep.) and the HERS CHEL Galactic plane survey HiGal 
(IMolinari et al. 2010 ) revealed a new population of very 
young sources, detectable at 70 /im, but yet dark at 24 



/im. Compa ring the embedded protostars found in Henning 



et al. (2010) to the starless clumps we found in Gil. 11, we 
conclude that one out of six starless clumps harbors an 
embedded source, that is invisible at 24 /im. This suggests 
that not all clumps presented here will be starless at 70 /im. 
Future studies will need to clarify the HERSCHEL view of 
our starless clumps. 



4. Distance independent results 

Using CLUMPFIND, we therefore extracted 901 clumps 
with peak column densities above our threshold of 1 x 
10 23 cm 2 . We found that 291 clumps have a Spitzer coun- 
terpart classified as YSO using the Gutermuth criteria. 
For 238 objects, STARFINDER identified a 24 /im (point) 
source in the MIPSGAL images, which had no YSO inside. 
During the visual inspection of the remaining 372 clumps, 
103 additional 24 /im sources were found, while 59 clumps 
were found to be partially or completely saturated in the 
MIPSGAL 24 /im images. Therefore, 210 clumps, or ~ 23%, 
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Sti 



.3 



15.0 14.0 
Galactic longitude / [ degrees ] 



Fig. 4. RGB image of the Galactic plane with Galactic latitude 1 = 10° to 20° using GLIMPSE 8 fim, MIPSGAL 24 



/j,m, and ATLASGAL 870 /im, respectively. Overplotted are CO contours from Dame et al. 
indicated as circles. 



(2001). Starless cores are 



show no signs of a heating source with the data employed. 
These can be considered as starless. 

The large number of visually identified 24 fim sources 
show that this step was crucial for a reliable source cata- 
log. Unfortunately, 24 fj,m point sources are often hidden 
in the unsteady background emission, hence algorithmic 
point source extraction is unable to distinguish the weak- 
est sources. The positions of the starless clumps are shown 
in a three color image, Fig. [4j and full details are given in 
Table [3j These clumps build a sample of potential starless 
clumps. As we show, most of them can be considered as 
massive. 



4.1. Results based on the classification 

It can be seen from Fig. [T] in the online appendix that al- 
most all of the clump s are embedded in larger structures (cf. 
Schuller et al.||2009| . These form filaments with lower den- 
sity gas, indicated by the 3 er contour. Only a few clumps 
seem to be isolated. 



The majority of the gas is concentrated towards the 
Galactic plane, with a small offset towards negative lati- 
tudes. Enhanced concentrations of gas/clumps are visible 
towards known regions, mainly W31, W33, M17, M16, and 
W39 (from west to east). 

A study to identify IRDCs, sol ely based on SPITZER 8 
fj,m extinction, was conducted by |Peretto fc Fuller| ([2009 ) . 
They found the fraction of starless IRDCs to be 32%. This 
is similar to the fraction of 23% found within this work. 



Nevertheless, as pointed out in Peretto & Fuller (20101, 
the detection of column densities via SPITZER 8 /Ltm ex- 
tinction breaks down at column densities larger than ~ 1 
x 10 23 cm~ 2 , which we require as minimum column den- 
sity in our study. In addition, extinction is very unlikely to 
be observed on the far side of the Galaxy. Therefore, more 
than a quarter of the complexes have no |Peretto &: Fuller] 
(2009) IRDC close by and are a completely new sample, 
which is likely at the far side of the Galaxy. 

The CLUMPFIND algorithm calculates the effective ra- 
dius of the clumps by equating the area of a theoretical cir- 
cular clump to the sum of the pixels. Results for the clump 
radii are shown in black in Fig. [5j The clump radii range 
between 10"and 40", with an average radius of 18". 




10 20 30 

effective angular radius / [ 



40 



Fig. 5. Histogram of the effective radius derived by 
CLUMPFIND for all clumps (black), for clumps with IRDC 
connection (hatched red area) , only, and for clumps without 
any IRDC indication (solid green area). 

4.2. Column densities 

The fluxes at the peak positions of the clumps can be used 
to derive a beam-averaged peak column density. Within 
their large NH3 survey, Wienen et al. (submitted) mea- 
sured the rotation temperature of 15 of our starless clumps 
directly. Both the mean and median temperature are T = 
15 K at these peak positions. Thi s is in agreemen t with 
temperature estimates for IRDCs (j Sridhara n et al. 2005 



Pillai et al.| [20061 |Vasyunina et ai.[ |2010[ |Peretto et "aT 
2010). Since we required all clumps to be devoid of 24 jim 



emission, we assumed that all our clumps have similar tem- 
peratures. We calculated the column density of the gas via 



Nh 2 = 



B\{\,T)m H2 KVL 



(1) 



for a gas-to-dust ratio of R = 100, where F\ is the flux at 
the given wavelength, B(A,T) the blackbody radiation as a 
function of wavelength and temperature, m# 2 the mass per 
H2 molecule, and fl the beam size. The mass absorption 
coefficient k = 0.77 cm 2 g -1 is adopted from the values 



given in Hildebrand ( 1983 ) using a dust emissivity index ft 
= 2 and an emissivity at 250 /zm of 3.75 x 10~ 4 . This is 
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Fig. 6. Logarithmic histogram plot of the column density 
(lower x-axis) and flux (upper axis). While the black his- 
togram represents the full sample of starless clumps, the 
hatched red and filled green histograms correspond to the 
near and far sample, respectively. (For details, cf. Sec.[5j) 



consistent with th e value for the diffuse ISM in Osscnkopf 
& Henning ( |1994[ ), a frequently used value for dark clouds. 
The calculated column densities for the starless clumps are 
given in table [3] For intermediate volume densities of 10 5 
cm and th in ice mantles, one can extra polate k = 1.85 
cm 2 g -1 from|Ossenkopf fc Henning ( 1994), as, e.g., used by 
Schuller et al.| ( |2009p . Including their different assumption 
of the mean molecular weight for t he ISM, column density 
estimates in | Schuller et al. ( |2009| ) would be smaller by a 
factor of three. 

The peak column densities vary only slightly. As shown 
in Fig. [6j 94% of the starless clumps have column densi- 
ties in the range between our survey threshold 1 x 10 23 
cm~ 2 and 2 x 10 23 cm~ 2 . Only 5 clumps, or 2%, have 
a peak column density larger than 3 x 10 23 cm~ 2 , which 
corresponds to 1 g cm" 2 . Nevertheless, the beam at a dis- 
tance of 3 kpc corresponds to 0.26 pc, hence is too large to 
resolve individual cores (~ 0.01 pc - 0.1 pc). The given col- 
umn densities are beam-averaged over large spatial scales 
and the actual peak column densities could be considerably 
larger. This effect preferentially reduces the column densi- 
ties of clumps further away more significantly than those of 
nearby clumps, introducing an artificial difference between 
the clumps on this side and on the far side of the Galaxy. 
This difference is clearly illustrated by the red and green 
histograms in Fig. [6] for the near and far clumps, respec- 
tively (see Sec. [HJ). 

To get a feeling for t he sm all-scale peak column den- 



sities, |Vasyunina et al.| (20091 studied the effects of dis- 
tance and telescope resolution onto the peak column den- 
sity. First, they produced an artificial r _1 density distribu- 
tion grid of 2000 AU, or 0.01 pc, resolution. Secondly, they 
smoothed the grid with different Gaussian kernels to im- 
itate observations with a 24" beam at different distances. 
They next compared the obtained column densities to the 
unsmoothed peak column densities and calculated from 
these ratios correction factors for different distances, which 




Fig. 7. Artist impression of face-on view of the Milky Way 
by R. Hurt (SSC-Caltech) / MPIA graphic. Plotted on top 
are the starless clumps presented here with the distance 
according to the distance flag in Table [3] 



resemble the peak column density seen with a linear spatial 
resolution of 2000 AU. The correction factors applied for a 
distance of ~ 2 kpc start at around ~ 17, and go up to ~ 
40 for distances around 4.5 kpc. Even assuming a minimal 
correction factor of ~ 10 and applying it to our sample, all 
clumps should contain smaller subregions of higher column 
densities, larger than 3xl0 23 cm~ 2 , following the Krumholz 
criterion for high-mass star formation. Nevertheless, this 
procedure cannot be applied to clumps at all distances. For 
clumps at the far side of the Galaxy in particular, the beam 
averages over several/many clumps and projection effects 
become more likely. 



5. Distances 

To determine additional physical parameters, in particu- 
lar the size and the mass of the clumps, the distance is a 
major parameter. A Galactic rotation curve was utilized 
to determine distances from the clumps' radial velocities. 
In the following, only clumps with a distance estimate are 
discussed. 

As the idea of this study was an unbiased survey of a 
large area of the sky with continuum data, a priori we have 
no information about the distances to the clumps found. 
To tackle this pro blem, we employed the Galactic rota- 
tion curve given in Reid et al. (2009). The necessary ve- 
locities are provided by Wienen et al. (submitted) . Wienen 
et al. (submitted) conducted spectroscopic follow-up obser- 
vations of NH3 towards bright peaks in the ATLASGAL 
survey. If no counterpart was found in Wienen et al. (sub- 
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Property Reviewed 


Near 


Far 


Only Far 


Estimated Error Near / Far / Only Far 


number of clumps 


115 


26 


9 




distance / [kpc] 


3.1 


12.8 


16.9 


0.5 


average effective radius / [pc] 


0.3 


1.0 


2.2 


0.28/1.2/1.6 


average mass / [M©] 


620 


5560 


26400 


factor of 4 


median mass / [M©] 


320 


4600 


21400 


factor of 4 


particle density / [cm -3 ] 


l.Ox 10 5 


2.3x 10 4 


1.3x 10 4 


factor of 2 



Table 1. Overview of typical clump properties for near and far clumps. The origin of the uncertainties is explained in 
Appendix [A] 



mitted), we used the HCO + survey of BGPS sources by 
Schlingman et al. (2011). 



To maximize the number of clumps to which we could 
assign a velocity, all clumps that lie within the same low- 
est significant contour were assumed to be connected. With 
this assumption, we were able to assign to connected clumps 
the same velocities as their neighbors. Incorporating all in- 
formation at hand, the velocities of 150 starless clumps, or 
71%, are known. The uncertainties in the velocities can be 
estimated by comparing the NH3 and HCO + velocities of 
clumps that have both measurements. We note that 134 of 
all clumps (not only starless clumps presente d here) were 



observed by both Wienen et al (submitted) and | Schlingman 
et all (2011 ). The average difference between both velocity 



measurements is 0.5 km/s, while their median difference 
is 0.3 km/s, with 2.3 km/s being the largest difference. 
Therefore, we estimated the velocity uncertainties to be 
0.5 km/s. 

Owing to the rotational structure and symmetry of the 
Galaxy, a Galactic rotation curve usually yields two dis- 
tance solutions for a given direction and velocity. To solve 
this distance ambiguity, additional information or assump- 
tions are required. While in studies of IRDCs it has often 
been argued that all sources lie at the near solution because 
their identification requires a bright mid-IR background, 
this argument could not be adopted here. The optically 
thin dust emission at 870 /im instead allowed us to identify 
clumps across the entire Galaxy. Nevertheless, coincidence 
with an IRDC favors the near solution and we used the cat- 



alog of IRDCs given in |Peretto fc Fuller| p009| to identify 
nearby objects within our sample. Although they co ver a 
different column density range (for details see section 4.1 ), 
considerable overlap can still be expected. During the visual 
inspection of the 24 /im emission, additional dark patches 
connected to our clumps were identified and noted as IRDC. 
In the following, all clumps with an associated IRDC were 
assumed to be on the near side. 

For 9 sources with velocities between -5 km/s and -1 
km/s, the rotation curve only allows the far solution, be- 
cause their near solution is meaningless (it places the source 
in the outer Galaxy, while we looked in the opposite di- 
rection towards the inner Galaxy). The far solution places 
them in, or clo se to, the Norma arm at ~ 17 k pc. However 
as disc ussed in Dame & Thaddeus ( |2008 ) and Green et al 



(2011), the velocities could also place them in the near 3 



kpc arm at ~ 5.2 kpc distance. For consistency with the 
adopted Galactic rotation model, we prefer the Norma so- 
lution. Future studies of these clumps could also use Hi self- 
absorption orJ^CO associations with well-known regions 



(Liszt et al. 1981) to better solve the distance ambiguity. 



clu mps a re on the far side of the Galaxy (cf. Table [T] and 
Sec. 6.2). Only few starless clumps on the far side of the 
Galaxy have been known previously (?), thus about a quar- 
ter of the sources are newly identified. Fig.[7]shows the loca- 
tions of the starless clumps within the Milky Way Galaxy. 
One notes a clear gap between 5 kpc and 11 kpc in the 
source distribution, which can be explained in several ways: 
(1) The elliptical orbits in the bulge of the Milky Way ran- 
domize its clouds' velocities and the rotation curve places 
them at random distances. (2) Circular orbits close to the 

tangent point have very large , hence small errors 

in the velocities propagate into large distance offsets. (3) 
The majority of the cold gas is homogeneously distributed 
in a molecular ri ng around the Galactic ce nter with 4 kpc 



< R G c < 8 kpc (Solomon & Rivolo 1989). Therefore, no 



clumps are expected outside that region 

6. Masses and clump mass function 

6.1. Masses 

Assuming optically thin emission, the mass of these clumps 
can be calculated from the dust continuum emission via 



ji/f " 

gas ~ B x (\,T)k' 



(2) 



where most of the parameters are the same as defined in 
Sec. |4.2| and d is the distance. Therefore, the mass can 
only be calculated for sources with distance measurements. 
For completeness, the mass is calculated for both the near 
and far solutions produced by the distance ambiguity and 
listed in Table [3j Fig. [8] shows the calculated masses for 
both the near and far solutions. The solid line indicates the 
theoretical sensitivity limit of our source extraction. The 
uncertainties in the masses are discussed in Appendix |A") 

We expect the amount of dense gas per volume to be 
similar on both sides of the Galaxy. We can therefore con- 
duct a consistency check by comparing the mass of the near 
and far population relative to the volume covered. 

All clumps on the near side of the Galaxy have masses 
below 10 4 Mq (cf. Fig.[8j black dots), while the maximum 
mass within the far-clumps is a factor of about three higher 
than the most massive clump on the near side. In addition, 
when looking at Fig. [8j the total number of near objects is 
clearly higher than the number of far objects. To make a 
quantitative comparison of both populations, we estimated 
the common sensitivity limit for distances between 10 kpc 
and 15 kpc to be 1 x 10 3 M Q . We next calculated the 
volumes that we cover on the nea r and the far side usin g 
the scale height of ^40 pc given in Bronfman et al. (2000) 



In summary, out of the 160 sources with velocity mea- 
surements 115 clumps are likely on the near side and 35 



Adding up all masses above the far detection threshold for 
the near range, kpc to 5 kpc, the total mass of clumps 
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distance / [kpc] distance / [ kpc ] 



Fig. 8. The mass in solar masses is plotted over the distance 
in kpc. The solid line indicates the sensitivity/completeness 
limit of our clump extraction, which depends on the dis- 
tance. The black filled dots represent sources with IRDC 
associations, and, therefore, the near distance is plotted. 
The red triangles represent the far solution of sources with- 
out an IRDC for which both a far and a near solution 
can be calculated; in addition, the green asterisks repre- 
sent the corresponding near solutions. Blue diamonds show 
sources with such low velocities that only a single solution 
can be found in the given direction. The error bar in the 
top right corner indicates uncertainties of a factor of two 
for the masses and 0.5 kpc for the distance. 



on the near side is ^37000 M Q . Doing the same for the far 
range 10 kpc to 15 kpc and normalizing it to the volume 
covered on the near side the total mass becomes 38000 M . 
Both masses agree given this crude estimate, supporting the 
idea that, statistically, the allocation is reliable. 

6.2. Observational biases for clumps on the far side 

The change in resolution over the survey's depth affects the 
sensitivity and the recognition of substructure significantly. 
The 19.2"beam corresponds to 0.28 pc linear spatial reso- 
lution at a distance of 3 kpc, and in contrast corresponds 
to 1.4 pc at a distance of 15 kpc. While the distances vary 
by more than an order of magnitude, the angular sizes of 
the extracted clumps vary by only a factor 2-3 and show 
no correlation with distance (cf Fig. [5| . This results in an 
almost linear correlation between the physical size and the 
distance, which is shown in Fig. [9} In addition, a single, un- 
resolved source would be 25 times fainter at 15 kpc than at 
3 kpc. As indicated by the solid line in Fig. [8j the complete- 
ness limit changes with distance. Both effects are studied 
in detail in the following. 

Taking an ATLASGAL map of 3°x2°, we re-extracted 
all clumps with CLUMPFIND using the same thresholds 
as explained above and calculated their masses assuming 
a generic distance of 3 kpc. In addition, we convolved the 
same map with a Gaussian profile, emulating a resolution of 
96", reducing the resolution by a factor of 5. This resembles 
the appearance of the same structure as seen at 15 kpc. 



Fig. 9. Effective radius in pc over the distance. Colors and 
symbols are as in Fig. [8] 



Once again CLUMPFIND was used to search for clumps 
using the same parameters, but assuming a distance of 15 
kpc when determining the mass. 

While 90% of the total mass was recovered in the lower 
resolution maps, the number of clumps extracted differed 
significantly. In the full resolution map, 328 clumps were 
extracted, whereas in the lower resolution map, only 20 
clumps were found. This implies that structures, which can 
be resolved into several clumps on the near side, cannot be 
resolved on the far side and that their fluxes then add up. 
The volumes covered at the two distances differ by a factor 
of about five. This would add to the probability of chance 
alignment and since the dust emission is optically thin, sev- 
eral faint clumps within the same beam may add up and 
may detectable. 

Therefore, one should keep in mind that clumps discov- 
ered on the far side are slightly different types of objects. 

6.3. Clump-mass function (CMF) 

Mass distributions connected to stellar populations or star 
formation are very often c ompared t o a p ower-law distri- 
bution, as first stated by Salpeter ( 19551 and discussed 
since then. The more recent study of |Reid et al. (20101 
even found that statistical errors allow the slope of (m) any 
sort(s) of astrophysical random distribution(s) to be fitted 
by a power-law with spectral indices similar to the Salpeter 
value of 2.35 within the errors. Nevertheless, a clump-mass 
function (CMF) is of great interest in connecting clumps 
to core- and star-formation efficiencies. To ensure that our 
data has sufficiently high quality statistics, we considered 
all 125 near clumps for a 'local' CMF and the 27 far clumps 
for a CMF on the other side of the Galaxy. As shown in Fig. 
[8j the completeness limit for the near and far populations 
is lower than 100 M Q and 1000 M@, respectively. The black 
histogram in Fig. [10] represents the number distribution of 
the near clumps dN/dM. Owing to the completeness limits, 
we cannot rely on the low-mass end of the CMF below 100 
Mq, nor do we have sufficient statistics to determine the 
existence of a broken power-law distribution. We instead 
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M /[ m ] 

Fig. 10. Clump mass function with the near population 
plotted in black, the far population plotted in red. Both 
populations' high-mass tail have been fitted, with the fitting 
range indicated by the thick horizontal bars at the bottom. 
The power law indices are fitted to -2.2 and -2.6 for the 
near and far population, respectively. 



only consider the high-mass tail. The green line indicates a 
fit to the high-mass tail of the near-clump mass function, 
with a logarithmic slope of a = -2.2. The red histogram 
is the far-clump mass function, fitted by the blue line of 
slope a = -2.6. Neglecting the uncertainties on the data 
points, the uncertainties on the slopes calculated by the 
IDL routine 'UNFIT' are 0.2 and 0.4 for the near- and far- 
clump mass functions, respectively. The thick bars at the 
bottom indicate the fitting range. If we expand the fitting 
range of the near population and include lower mass bins, 
the slope becomes immediately shallower. If we reduce the 
fitting range, the slope remains constant within the uncer- 
tainties. Therefore, we believe that the fitting range is re- 
liable and the uncertainty is reasonable. For the far range, 
the situation is more difficult because of the smaller num- 
ber of bins. An enlargement of the fitting range to smaller 
bins again reduces the slope dramatically. Including the 
next mass bin does not change the result significantly, but 
including additional bins of lower mass would steepen the 
slope to -3.1. This number is almost within the errors and 
can be entirely explained by the lower quality statistics. 
Nevertheless, the fit to the far mass distribution should be 
interpreted with caution. 



7.1. Mass of the most massive star 

To place constraints on the lifetime of starless clumps, we 
first need to estimate what clump mass is required so that 
the final cluster can house at least one massive star. 

Star-formation efficiencies on scales from clumps to 
stars do not have a common v alue but a number of studies 
estimate that it is 23% - 50% ( |Chabrier fc Hennebelle|2010 
and references therein). For high-mass star formation, the 
numbers are even more weakly constrain ed (5% - 50%, cf. 
Krumholz et~aT|2007[ |Kuiper et aL|2010[). 



Following the definition given in Williams et al. (2000) 
and Beuther et al. (2007), these clumps will most likely 
form entire clusters instead of single stars. To estimate the 
required clump mass to form a star of given mass, one must 
assume a gas-to-star formation efficiency and the initial 
mass function (IMF) of the cluster produced. 



We assume the IMF of Kroupa (2001 ) and normalize it 



to the probability that at least one star with a mass higher 
than 20 M is formed. Integrating the normalized IMF over 
the expected mass range of stars, from 0.08 to 150 M Q , the 
stellar mass of this cluster is on the order of 320 M . With 
a star formation efficiency (SFE) of 30% (for details see 
next paragraph) , we estimate the mass of a clump with the 
potential to form at least one star more massive than 20 
Mq to be 10 3 M Q . A 3 x 10 3 Mq clump is required to 
form a 40 M Q star (see also Table [2j) . 

Since the mass distribution follows a power law, the 
number of clumps with masses higher than a given thresh- 
old is very sensitive to that threshold. In the picture in 
which these clumps form entire clusters following the IMF, 
the estimate of the stellar cluster mass relative to the most 
massive star seems quite reliable. Nevertheless, the esti- 



mates of the SFE vary over a wide rang e (|Lada fc Lada 
20031 [Alves et al.[ [20071 |Johnston et~aLl |2009[ |BontelrIps 



et aH[2010bD. Here, we us ed th e SFE given in Lada fc Lada 
( 2003), | Alves et al.| ( |2007| ) , and |Bontemps et al.| ( |2010bf ot 
30%"^ 

In the near sample derived across 20 deg 2 of the sky, 
with this estimate only 14 starless clumps have the poten- 
tial to form stars more massive than 20 M 9 , and only 3 
have the potential to form a 40 M Q star. 

The MIPSGAL images of the 3 most massive clumps 
are shown in Fig. 11 As can be seen, none of the regions 



are isolated but all are connected to regions already con- 
taining 24 /jm sources. At the distances to these objects 
varying between 3.6 kpc and 4.6 kpc, their effective radii 
become between 0.7 pc and 0.8 pc. Therefore, their particle 
densities are not among the highest of our sample at only 
4 x 10 4 cm~ 3 to 8 x 10 4 cm~ 3 . Nevertheless, these are 
very interesting objects and very promising in the context 
of massive star formation. 



7. Lifetimes 

As discussed in Sec. |6.21 clumps seen at the far side of the 
Galaxy are a mixture of clumps seen in chance- alignment 
or groupings that we would identify as several structures on 
the near side. To form a more consistent sample, in the fol- 
lowing we only consider starless clumps that are identified 
on the near side. 



7.2. Lifetime of clumps 

Since there is no reason to believe that the starless clumps 
have ages that are correlated with those of other objects in 
this sample, the sample should span the entire age range 
expected for these clumps. If this is the case, the oldest 
clumps will start forming stars as new and similar clumps 
appear. Therefore, their lifetime can be calculated by com- 
paring the number of these clumps identified to the number 
of massive stars formed. 
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SFR 

Mass of Potential Star 


H83, diffuse 
20 M 


ISM opacity 
40 M Q 


OH94, cold d 
20 M 


ense opacities 
40 M 


Clump-mass Threshold 
Number of Clumps Above Threshold 


1065 M 
14 


2960 M 
3 


1065 M 
6 


2960 M 
1 


SFR of 1 M / yr 
SFR of 3 M Q / yr 
SFR of 6 M Q / yr 


3.3 x 10 b yr 
1.1 x 10 5 yr 
5.5 x 10 4 yr 


2.0 x W b yr 
6.6 x 10 4 yr 
3.3 x 10 4 yr 


1.4 x 10 b yr 
4.8 x 10 4 yr 
2.4 x 10 4 yr 


6.6 x 10 4 yr 
2.2 x 10 4 yr 
1.1 x 10 4 yr 



Table 2. Lifetimes of starless clumps calculated for different sets of parame ters. Estimates are calculated using different 
Milky Way star formation efficiencie s for opacities for both the cold ISM (Hildebrand 1983| H83) and dense but cold 
regions ( Ossenkopf fc Henning||1994| OH94). 




Fig. 11. MIPSGAL images of the three most massive starless clumps found in this survey, with (from left to right ) 7400 
Mg, 4700 Mq, and 3000 Mg. Green contours are from ATLASGAL, red lines mark the boarders of the clumps identified 
by CLUMPFIND. 



To do so, we need to 1) extrapolate the number of 
starless clumps we would find in the entire Milky Way 
galaxy, and 2) estimate the number/fraction of massive 
stars formed every year: 

1) We assume that most star-forming gas is distributed in 
a ring around the Galacti c center, between 4 kpc a nd 8 kpc 
Galactocentric distance ( Solomon fe Rivoiol 19891 with a 
scale height of 40 pc (Bronfman et al. 2000). Taking into 
account that we only consider clumps up to a distance of 
5 kpc from the Sun within the direction of the survey, and 
above the clump mass thresholds for 20 M and 40 M 
stars, for the whole Milky Way Galaxy, we expect to iden- 
tify 1043 and 223 clumps, respectively. 

2) To calculate the fraction of massive stars formed ev- 
ery year, one needs to assume a star-formation rate (SFR) 
for the Milky Way. The most recent publica tions suggest 
a star formatio n rate of around 1-2 M Q / yr ( |Robitaille fc 
Whitney[ |2010). 



For a SFR of 1 M Q /yr, the lifetimes of massive starless 
clumps become 3 x 10 5 yr and 2 x 10 5 yr calculated for 
all clumps more massive than 10 3 M Q , and 3 x 10 3 Mg, 
respectively. For higher SFRs, the lifetimes become shorter 
in a linear fashion, e.g. for 6 M /yr the lifetimes become 
5.5 x 10 4 yr and 3.3 x 10 4 yr, respectively. Lifetimes for 
different parameters are summarized in Table [2j As we dis- 
cuss in Sec. \8.2\ we estimate the lifetime of these objects to 
be (6 ± 5) x 10 4 yr. 



8. Discussion 

8.1. Discussion of the clump-mass function 

The slope of the CMF for starless clumps found in this 
study, a — -2.2, is similar to the value of the Salpeter IMF. 
Nonetheless, as we study clumps that will most likely host 
small clusters rather than individual stars, we do not pro- 
pose a one-to-one mapping to the IMF. We emphasize that 
not all clumps will eventually form clusters or even be tran- 
sie nt objects. Another set o f massive clumps was presented 



in |Peretto fc Fuller ( 2009 ) , including a subset of clumps 
without MIPS 24 /im emission. They found the mass func- 
tion of IRDCs to be simi l ar to the CO clump mas s distribu- 
tion (|Simon et al.[ [20061 IPeretto & Fuller! |2010|) with a 
-1.7. In 



addition, they used a derivative of CLUMPFIND 
to search their extinction maps for structures within the 
IRDCs. These fragments have a rather Sa l peter- like slope 



2006 Peretto & Fuller 



2010) which is 



Although the extinction method 



( Rathborne et al. 
then si milar to our result 

used by Peretto & Fuller ([2009 ) to calculate column densi- 
ties of IRDCs is sensitive to a lower column density range 
than that of the starless clumps we present, the 870 /^m 
emission iden tifies objects similar to fragments |Peretto ~&\ 
Fuller] ([2009| find in IRDCs. 



Th e studies 
p005], and 



of 



Williams et al.| 
(2006; 



Reid fc Wilson 
all targeted the more 



Beltran et al.| 
evolved high-mass protosteliar objects. They found a com- 
mon break in the CMF at 100 M and fit power laws to 
their high-mass end with exponents between -2.0 and -2.32. 
Although we did not attempt to fit the exact break, a break 
point of 100 M Q or just above seems to be in good agree- 
ment with the starless CMF, but that clearly might be bi- 
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ased by the completeness limit close to 100 M Q . However, 
the exponent of the near CMF agrees with all values within 
the errors. Thus, comparing our results to earlier studies of 
more evolved clumps shows that there is no evidence that 
the CMF of starless clumps is different from a CMF at later 
evolutionary stages. 

Thus, we found that the exponent of the CMF for 
clumps at the far side of the Galaxy is larger than most 
other values. Although the uncertainty is larger and the dif- 
ference could be explained by the errors, this trend might 
equally repre sent a general scale-dependent trend. |Beltran| 
et al. ( |2006[ ) distinguished between the populations at d 
< 2 kpc and d < 6 kpc and found that the exponent for 
the more clearly resolved population < 2 kpc is shallower 
than for d < 6 kpc. Our far population is even more distant 
than their sample and the far CMF's exponent would con- 
tinue the trend to steeper slopes. This could be interpreted 
as a general scale-dependent trend and rather a matter of 
resolution than of true structure. 

Studies of the core mass function for low-mass star for- 
mation target significantly smaller objects. To have suf- 
ficient spatial resolution, they are typically chosen to be 
nearby. Nevertheless, whether they present starless cores or 
more evolved objects, many studies have measured power- 



1998 



law slopes close to the Salpeter value ( Motte et al. 
Johnstone et~al| [20001 |Alves et all |2007[ |Enoch et al. 
2008). 



8.2. Discussion of the lifetimes 

The estimates of the lifetimes of starless clumps are based 
on the SFR of a given mass range and the number of clumps 
above a corresponding clump mass. Both estimates involve 
several assumptions, which may introduce errors. 

Lifetimes are inversely p roportional to the SFR. The 
SFR varies from 1 M^/yr dRobitaille fc Whitney) |2010[) 



to over 4 M Q /yr ( Diehl et al. " 2006p to even larger values 
but the most recent publications favor the smaller values 
Nevertheless, our own survey looks at a region in the vicin- 
ity of the Galactic center for which both |Rosolowsky et al. 
(2010) and ? found that more dense gas is located within 
the inner 1 < 30° of the Galactic plane. While most gas is 
concentrated within a molecula r ring around the Galacti c 



center at 4 kpc < R < 8 kpc (Solomon & Rivolo 1989) 



the ring does not seem homogeneous and the outer regions 
of that ring seem to contain less gas. This might indicate 
that the extrapolated numbers for the entire Galaxy might 
be higher than average and therefore require a higher SFR. 
This would reduce the lifetime estimate. 

The lifetimes are also proportional to the number of 
clumps of given m ass and when we use the ex trapolated 
dust opacities from Ossenkopf fc Henning| ( |1994[ ) our clump 
masses decrease and the number of clumps above our 
threshold is lower. This, in turn, reduces the lifetimes. In 
addition, one should keep in mind that small number statis- 
tics are involved. If one of these objects turns out to already 
be star-forming, the estimated lifetimes are reduced, while 
correcting chance alignments of mid-IR sources has the op- 
posite effect. 

Another factor contributing to the uncertainties are 
the star-formation efficiencies, which directly influence the 
mass thresholds themselves. As explained in Sec. |7.1[ the 
number distribution decays as a power of the mass, which 
makes the lifetime estimates very sensitive to the clump 



mass thresholds. On the basis of the range of possible life- 
times shown in Table [2j we estimated the uncertainties to 
be one order of magnitude. 

In summary, most effects seem to reduce the lifetimes. 
Based on these arguments, we estimate the lifetime of star- 
less clumps to be on the order of (6 ± 5) x 10 4 yr. 

For volume-averaged particle densities of 10 5 cm~ 3 , the 
free-fall time becomes ~ 1.6 X 10 5 yr, and lower densities 
increase this number. Therefore, the free-fall time is about 
a factor of two longer than the lifetime we found for starless 
clumps, but both agree within the errors. 



In good agreement with previous studie s (Motte etal 

2007] |Hatchell fc Fuller) [2008} |Motte et al.||2010| ), Table |2] 
shows that the estimated lifetimes of more massive clumps 



clumps 

are smaller. Accordingly, it was expected that the lifetimes 
we found for starl ess dense clumps are shorter than the 3 
x 10 5 yr found by Kirk et al. ( 2005 1 for low- mass cores. 



An adequate comparison to the lifetime of high-mass 
starless clumps is difficult because only a few studies exist. 
Most deal with s l ightly different kinds of objects. For exam- 
ple, [Motteetar ( 2007[ ) studied the nearby Cygnus X region 
and found cores that cover the same volume densities as our 
clumps, while thei r masses are significantly lower (see also 
Motte et al.||2010 ). However, they found clumps of similar 
mass to our objects but with lower volume densities. Using 
SiO as an additional tracer of star formation, IMotte et al.l 



(2007) did not find a single starless massive dense core in 
their sample. Therefore they proposed that the low-density 
starless clumps dynamically evolve into star-forming mas- 
sive dense cores. Their lifetime estimate of starless massive 
dense cores becomes < 10 3 yr. 

A similar study covering more simil ar entities both in 
density and mass by Russeil et al. ( 2010[ ) found one starless 
clump, hence they derived a statistical lifetime of ~ 1 x 10 4 
yr. This agrees with our estimate, but at the lower edge 
of the error. The difference could also be caused by their 
inclusion of SiO as a tracer of star formation, which could 
reduce the number of starless clumps we found. 

8.3. Comparison to other surveys 

8.3.1. Comparison to the Bolocam Galactic Plane Survey 
(BGPS) 

The Bolocam Galactic Plane Survey (BGPS) performed a 
systematic study of the inner Galactic plane at 1.1 mm. 
Although their rms across the region 10° < 1 < 20° is ~ 25 
mJy/beam, their overall 5 a level is 0.4 Jy, at which they 
have a completeness level of 99%. Their 0.4 Jy complete- 
ness threshold converts into a column density threshold of 
~ 6 x 10 22 cm -2 over their 31"beam when using the same 
assumptions as for our data, while their sensitivity trans- 
lates into a rms of 64 mJy/beam at 870 /im and their beam 
width of 31". With this column density, they found 1211 
sources on the same region as we studied here. 

For that data set, Dunham et al. (2011) used different 
point source catalogs to search within all BGPS sources 
for mid-infrared tracers of star formation activ ity. As re- 
sources they employed the Red Sources (R08, Robitaillc 



et al. 2008), the RMS catalog, the catalog of extended 



green objects (EGOs, [Cyganowski et al.|2008 ), and the full 



GLIMPSE catalog. The EGOs have no direct counterpart 
in this study, but trace shocked gas. To be truly starless, 
starless clumps should not contain EGOs. 
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Reviewing the Dunham et al. (20111 source catalog for 
the same region as surveyed in this study, we found that 
for ~ 70 % of the sources they found neither RMS sources, 
R08 sources, nor EGOs. 



Following a similar approach to Dunham et al. (2011 1 
and using the identical mid-IR tracers, the R08, RMS, 
and EGO catalogs, for the ATLASGAL clumps, we found 
R08 counterparts within 163 clumps, 39 RMS counterparts 
within clumps without R08 sources, and 9 matches among 
EGOs and clumps without the previous tracers. This leaves 
719, or 77%, of our ATLASGAL clumps starless. Both frac- 
tions of infrared quiet clumps a re significantly larger than 
those of |Peretto fc Fuller | ( [2009] ) (32%) or we found (23%) 
including all tracers. Even without visual inspection but in- 
cluding MIPS GAL 24 urn sources found by STARFINDER, 
we would have found only 42% of the clumps to be starless. 
These differences in our statistics from both of the latter 
studies can be explained by the superior sensitivity of the 
"by eye" source confirmation and the higher spatial resolu- 
tion of the MIPSGAL survey over the MSX images. 



8.3.2. Comparison to PLANCK's Early Cold Cores 

That Planck's Early Cold Core Source List (ECC) (?) does 
not contain any cold sources in the region of the sky we 
surveyed here is apparently remarkable. However, the ECC 
contains only sources colder than T < 14 K averaged over 
PLANCK's beam sizes of several arc minutes. Owing to 
the high gas density and ongoing star formation as well as 
confusion, no sources are expected to be found within the 
Galactic plane (?). 



9. Conclusion and outlook 

9.1. Conclusion 

For the first time, we have presented a complete and un- 
biased sample of high-mass starless clumps on 20 deg 2 of 
the sky. To concentrate on the actual potential precursors 
of massive stars, we imposed a minimum peak column den- 
sity of 1 x 10 23 cm~ 2 . From ATLASGAL, we extracted 
901 clumps across the region 10° < 1 < 20° of the Galactic 
plane. Using the GLIMPSE source catalog and MIPSGAL 
24 fxxsx images, we compared clumps found at 870 /jm to 
near- and mid-infrared tracers of ongoing star formation. 

Neglecting clumps that are saturated at MIPSGAL 24 
fim, we identified 210 clumps, or 23%, to be starless. Their 
effective radii range from 10" to 40", and most of the beam- 
averaged peak column densities are in the range 1 - 2 x 10 23 
cm~ 2 . Correcting the single-dish peak column densities to 



true p eak column densities as discussed in |Vasyunina et aL 
( 2009 ) even suggests that all clumps should have peak col- 



umn densities above the fragmentation thres hold for mas- 
sive star formation of 1 cm proposed by Krumholz & 



McKee| ( |2008 ). 

For ~ 71% of the starless clumps, we had the veloci- 
ties and were able to calculate their distance. If a clump 
is connected to an IRDC, then we assumed a distance on 
the near side. We found that about a quarter of the starless 
clumps lie on the far side of the Galaxy and were therefore 
previously unknown. 

We found that the masses of starless clumps on the near 
side range from 10 M to 7500 M , and that objects on the 
far side have masses between 1000 M Q to several 10000 M . 



The different mass regimes are a consequence of our lim- 
ited spatial resolution, biasing us towards detecting larger 
structures as a single clump on the far side, which we would 
resolve into several clumps on the near side. The mass dis- 
tributions of the near and far populations could be fitted by 
power laws with slopes a = -2.2 and a — -2.6, respectively, 
and agree within the errors. This shows that the mass dis- 
tribution of clump populations on the near and far side of 
the Galaxy do not differ significantly. 

Drawn from the population of clumps on the near side, 
we found that only 14 objects are massive enough to form 
clusters with stars more massive than 20 M . Only 3 star- 
less clumps have the potential to form stars even more mas- 
sive than 40 M . We estimate the minimum clump mass 
required to form a cluster with a 20 M or 40 M star 
to be 10 3 M or 3 x 10 3 M , respectively. Since the star- 
formation efficiency used is an upper limit, these numbers 
are lower limits. 

Extrapolating the numbers of massive starless clumps 
from our survey volume to the Milky Way Galaxy, we es- 
timated the lifetime of the most massive starless clumps 
to be on the order of (6 ± 5) x 10 4 yr. To do so, we as- 
sumed a star formation rate of 1 M to 3 M per year for 
the entire Galaxy. We also discussed a possibly enhanced 
star formation activity within the surveyed volume and its 
implications for the assumed global star formation rate. 

9.2. Outlook 



In the future, we wish to extend our classification of the 
clumps to more evolved stages. To identify massive young 
stellar objects in particular, the Red MSX Source (R MS) 
survey will be used. Incorporating CORNISH (see Sec. 3.4) 
and RMS, in addition to GLIMPSE and MIPSGAL, we 



plan to establish an evolutionary sequence of the early 
stages of massive star formation. With this classification, 
we will be able to derive the relative timescales for the full 
sequence. Using the absolute timescale established in Sec. 
|7.2[ we will be able to translate the relative into absolute 
timescales and compare these to earlier studies. 

In addition, the objects that we have presented herein 
represent an ideal and unique sample for follow-up studies 
of the early stages of massive sta r formation. The upc om- 
ing HI-GAL /HERSCHEL survey ( |Molinari et aT|[20To| will 
help to provide additional constraints on the nature of these 
clumps, but deep integrations with the HERSCHEL 100 
fim band will be needed to uncover all embedded sources. 
In addition, the high resolution of ALMA will allow de- 
tailed studies of star-forming regions on the other side of 
the Galaxy. 
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Table 3. Properties of starless clumps as described in the paper. The full table is available in the online version of this paper. 
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Galaxy. However, their galactic latitude and velocities are also consistent with the near 3 kpc arm at a distance of 5.2 kpc (Dame & Thaddeus 
that Wienen et al. (submitted) have observed a position within the clump's boundary definition. If the NH3 fla g is absent but a JNH3' 
from neighboring clumps (see Sec.[H]for details). The same yields for the HCO + flag and HCO + observations in Schlingman et al 
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Appendix A: Errors and uncertainties 

As conventional error propagation breaks down when the 
uncertainties become larger than a few percent, one can 
only point out the individual sources of errors and estimate 
the final uncertainties. 

The distance error mainly stems from uncertainties in 
the Galactic rotation curve and errors in the gas velocities 
can be neglected. Error propagation inclu ding the velocity 
uncertainties following Reid et al. (2009) suggests uncer- 
tainties of smaller than 0.1 kpc. However, owing to intrin- 
sic errors and deviations from the global Galactic rotation 
we estimate the distance to be uncertain to within 0.5 kpc. 
This leads to a contribution to the final mass uncertainties 
ranging from 10% to 50%, that depends on the absolute dis- 
tance. In addition, individual objects close to the Galactic 
center that have non-circular orbits may be placed at ran- 
dom distances and contaminate the sample. 

In the literature, temperatures of sta rless cores range 



from 10 K to 20 K with the bulk at 15 K ( Sridharan et al. 
Pillai et al. 2006| |Vasyunina et a 



2005 



et al. 



20101 IPeretto 



2010). The temperature estimate here is based on di- 



rect observations of 15 out of 210 starless clumps and is in 
good agreement with earlier studies. A temperature uncer- 
tainty of ± 5 K at 15 K may introduce mass uncertainties 
of about a factor of two. 

The dust properties and the gas-to-dust ratio are very 
uncertain as well and might contribute another factor of 
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two to the errors. The flux uncertainties are domi nated by 



the calibration uncertainties, which are ~ 15% (Schuller 



et al. J2009|! 
hen 



When calculating the column density as well as the 
masses, the predominant uncertainties are those of the dust 
properties and temperatures. For the mass, the uncertainty 
in the distance is equally important. Altogether, the total 
uncertainties in the mass may be as large as a factor of 
five. 



Appendix B: Omitted regions 
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Table B.l. Regions where MIPSGAL 24 /zm images are to saturated to do a classification. Regions listed below were 
omitted during the classification. 



global identifier 


Gal Lon 


Gal Lat 


Ra 


Dec 


Peak Flux 


Total 
Flux 


Radius 




[°] 


[°] 


[°] 


[°] 


[ Jy / beam ] 


[Jy] 


["] 


8 


10.3225 


-0.1608 


272.2570 


-20.0861 


4.40 


45.36 


50. 


45 


10.1325 


-0.3775 


272.3610 


-20.3573 


1.30 


6.38 


28. 


59 


10.1392 


-0.3658 


272.3535 


-20.3458 


1.06 


5.59 


26. 


65 


10.1458 


-0.3158 


272.3103 


-20.3157 


1.00 


7.15 


31. 


70 


10.1808 


-0.3692 


272.3781 


-20.3109 


0.94 


6.48 


29. 


79 


10.1575 


-0.3775 


272.3739 


-20.3354 


0.88 


3.34 


23. 


96 


10.1442 


-0.3558 


272.3468 


-20.3365 


0.78 


2.34 


19. 


465 


14.9838 


-0.6956 


275.0954 


-16.2465 


6.60 


166.25 


74. 


467 


15.0121 


-0.7056 


275.1185 


-16.2262 


5.24 


118.69 


64. 


502 


14.9755 


-0.7139 


275.1082 


-16.2624 


1.90 


27.56 


49. 


514 


15.0022 


-0.7255 


275.1320 


-16.2444 


1.45 


17.96 


44. 


519 


14.9938 


-0.7305 


275.1325 


-16.2541 


1.34 


12.58 


36. 


554 


14.9738 


-0.7405 


275.1318 


-16.2765 


0.92 


4.37 


25. 


561 


14.9755 


-0.7372 


275.1296 


-16.2734 


0.91 


6.55 


28. 


617 


14.9688 


-0.7389 


275.1279 


-16.2801 


0.72 


2.42 


20. 


656 


13.9908 


-0.1208 


274.0768 


-16.8488 


0.61 


1.40 


16. 


687 


12.7975 


-0.2275 


273.5786 


-17.9485 


0.56 


2.69 


23. 


707 


13.9892 


-0.1358 


274.0897 


-16.8574 


0.54 


1.52 


18. 


786 


14.9738 


-0.7539 


275.1441 


-16.2827 


0.50 


1.33 


17. 


1114 


15.0545 


-0.6256 


275.0659 


-16.1512 


3.59 


28.11 


40. 


1116 


15.0511 


-0.6423 


275.0795 


-16.1620 


2.98 


46.09 


50. 


1119 


15.0678 


-0.6140 


275.0617 


-16.1339 


2.61 


22.05 


37. 


1122 


15.0995 


-0.6889 


275.1461 


-16.1414 


2.25 


22.50 


38. 


1125 


15.0778 


-0.6073 


275.0605 


-16.1220 


2.02 


21.50 


36. 


1126 


15.1844 


-0.6223 


275.1266 


-16.0351 


1.96 


37.36 


58. 


1131 


15.1311 


-0.6706 


275.1448 


-16.1048 


1.74 


21.13 


46. 


1137 


15.1028 


-0.6573 


275.1187 


-16.1235 


1.40 


31.14 


56. 


1139 


15.0028 


-0.7206 


275.1277 


-16.2415 


1.36 


9.88 


30. 


1140 


14.9929 


-0.7322 


275.1335 


-16.2558 


1.36 


11.98 


35. 


1141 


15.1061 


-0.6939 


275.1539 


-16.1378 


1.36 


12.79 


36. 


1152 


15.1111 


-0.7122 


275.1732 


-16.1421 


1.16 


11.03 


36. 


1153 


15.1161 


-0.6340 


275.1038 


-16.1008 


1.16 


9.14 


32. 


1157 


15.0961 


-0.7106 


275.1643 


-16.1545 


1.10 


5.29 


26. 


1168 


15.0944 


-0.7139 


275.1666 


-16.1576 


0.95 


3.61 


22. 


1174 


15.0944 


-0.6073 


275.0687 


-16.1073 


0.92 


4.57 


23. 


1180 


15.1078 


-0.6273 


275.0936 


-16.1050 


0.91 


5.14 


25. 


1181 


15.1061 


-0.5923 


275.0606 


-16.0900 


0.90 


5.60 


29. 


1187 


15.0944 


-0.7322 


275.1834 


-16.1662 


0.85 


4.65 


27. 


1193 


15.0412 


-0.6189 


275.0532 


-16.1598 


0.81 


5.39 


28. 


1202 


15.1278 


-0.6923 


275.1630 


-16.1180 


0.76 


3.21 


23. 


1207 


14.9845 


-0.7455 


275.1417 


-16.2694 


0.74 


2.95 


22. 


1214 


15.0995 


-0.7322 


275.1859 


-16.1618 


0.73 


2.98 


22. 


1215 


15.0778 


-0.6823 


275.1293 


-16.1573 


0.72 


4.79 


27. 


1223 


15.1028 


-0.6223 


275.0865 


-16.1070 


0.71 


3.04 


22. 


1224 


14.9845 


-0.6589 


275.0621 


-16.2286 


0.71 


4.20 


26. 


1231 


15.1278 


-0.7023 


275.1722 


-16.1227 


0.68 


1.34 


16. 


1233 


15.0861 


-0.6273 


275.0829 


-16.1241 


0.68 


2.48 


21. 


1285 


15.1228 


-0.7106 


275.1774 


-16.1310 


0.54 


1.19 


15. 


1291 


15.0561 


-0.6989 


275.1340 


-16.1843 


0.53 


1.55 


18. 


1294 


14.9879 


-0.7555 


275.1525 


-16.2712 


0.53 


0.87 


14. 


1297 


15.0328 


-0.7372 


275.1577 


-16.2229 


0.52 


1.78 


19. 


1301 


15.1094 


-0.5823 


275.0531 


-16.0823 


0.52 


1.33 


17. 


1302 


15.1244 


-0.7023 


275.1706 


-16.1256 


0.52 


1.18 


15. 


1303 


15.0245 


-0.7389 


275.1552 


-16.2310 


0.52 


0.68 


12. 


1315 


15.1311 


-0.6856 


275.1585 


-16.1119 


0.52 


2.16 


21. 


1324 


15.1011 


-0.6323 


275.0949 


-16.1132 


0.51 


1.64 


18. 


1327 


15.0628 


-0.6023 


275.0486 


-16.1328 


0.51 


1.82 


19. 


1331 


15.0278 


-0.7339 


275.1522 


-16.2257 


0.51 


1.68 


18. 


1737 


19.0742 


-0.2742 


276.6897 


-12.4356 


0.52 


2.73 


23. 



Notes: Columns are identifier, galactic longitude, galactic latitude, right ascension, declination, peak flux, integrated flux and 
radius as calculated by CLUMPFIND. 
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Appendix C: Stamps of starless regions 




Fig. C.l. MIPSGAL 24 image with ATLASGAL con- 
tours on top. Starless clumps are marked with a red aster- 
isk. The numbers correspond to the global identifier given 
in Table El 
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Fig. 1. ATLASGAL 3 and 6 a contours in yellow and green, respectively, on top of a 24 /im MIPSGAL image in 
logarithmic scale, with clumps plotted in overlay. Plus signs represent clumps for which the near solutions is assumed 
while circles are clumps with far solution assumed. For clumps plotted with a diamond only the far solution exists, while 
for sources with a box no velocity information is present. White stripes at the edges are artifacts from the MIPSGAL 
coverage. 
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Fig. 2. continued. 
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Fig. 3. continued. 
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Fig. 4. continued. 



